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INTRODUCTION 


This  addendum  details  the  statistical  procedures  utilized  In  prepar- 
ing the  report,  Decision  Risk  Analysis  for  XM204,  105MM  Howitzer,  Towed 
Reliability/Durability  Requirements  (PAA-TR-1-73) A computer  program 
was  developed  for  this  study  to  simulate  DT/OT  II  testing,  analysis  of  test 
data,  the  decision  made  based  on  this  analysis,  and  the  funding  implica- 
tions of  this  decision. 

A two-parameter  Weibull  family  was  assumed  to  describe  durability 
failures.  The  test  results  were  used  to  estimate  the  two  parameters 
for  each  subsystem:  carriage,  recoil  system,  tube,  and  breech.  This 

was  followed  by  a test  of  hypothesis  which  tested  whether  sufficient 
information  was  available  to  reject  the  hypothesis  that  the  required 
durability  had  been  obtained.  The  subsystem  durability  was  then  used 
to  compute  replacement  requirements  over  the  system  lifetime. 

PROBLEMS 

The  first  statistical  problem,  parameter  estimation,  was  unusual  in 
that  a large  number  of  subsystems  would  survive  testing.  This  occurred 
as  the  carriage  was  a long-lived  subsystem  and  the  test  truncation  point 
was  near  the  expected  carriage  life.  Also,  other  failed  subsystems  would 
be  replaced  until  one  of  the  following  occurred: 

a.  The  carriage  failed. 

b.  The  truncation  point  was  obtained. 

c.  No  additional  back-up  subsystems  were  available. 

Thus,  for  each  subsystem  the  test  data  would  consist  of  rounds-to-failure 

Hlazza,  Thomas  N.  and  Banash,  Robert  C. , PAA-TR-1-73.  Decision  Risk  Analysis 
for  XM204,  105MM  Howitzer,  Towed  Reliability/Durability  Requirements.  Systems 
Analysis  Division,  Plans  and  Analysis  Directorate,  US  Army  Weapons  Command, 
April  1973. 
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and  random  truncation  points.  Parameter  estimation  was  performed  by 
maximizing  the  log  likelihood  equations. 

The  next  problem  was  to  use  these  data  in  a test  of  hypothesis. 

This  was  accomplished  by  using  the  large-sample  distribution  properties 
of  maximum  likelihood  estimaters,  i.e.,  that  the  parameters  have  a joint 
normal  distribution.  The  a-percentile  confidence  interval  is  approxi- 
mately ellipsoidal;  the  a-percentile  durability  boundary  was  obtained 
by  maximizing  over  this  ellipse. 

The  third  problem  was  solving  for  the  number  of  renewals  in  a 

specified  time  interval.  This  was  accomplished,  initially,  by  the 

2 

method  of  Lomnicki  and  later  (to  save  computer  time),  by  taking  averages. 

The  purpose  of  this  addendum  is  to  detail  these  procedures. 

PARAMETER  ESTIMATION 

The  form  of  the  density  function  for  the  two-parameter  Weibull 
distribution  used  in  this  study  is: 

f(t)  = Aat  e dt  t ^0;a>0;A>0  . (1) 

Consider  the  following  subsystem-life  test.  M-like  systems  are  placed 

on  test  and  each  system  is  composed  on  one  critical  subsystem  and  several 

different  noncritical  subsystems.  A total  system  configuration  is  required 

to  conduct  the  test;  however,  with  respect  to  probability  of  failure, 

each  subsystem  is  assumed  to  be  independent.  As  each  non-critical  sub- 
system failure  occurs,  the  failure  time  is  noted,  and  the  failed  subsystem 

2 

Lomnicki,  Z.  A.,  A Note  on  the  Weibull  Process.  Biometrika  53,  375-381 
(1966). 
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is  replaced  with  an  identical  new  subsystem.  Each  of  the  M systems  con- 
tinue with  the  test,  with  failed  noncritical  subsystems  being  replaced, 
until  either  a predetermined  system  truncation  point  (for  example,  when 
a stated  number  of  rounds  have  been  fired)  is  reached  or  until  a critical 
subsystem  failure  occurs.  This  test  differs  from  the  well-known  Type  I 
and  Type  II  censoring  in  that  the  subsystem  truncation  points  are  ran- 
dom variables  as  Is  the  number  of  each  type  of  noncritical  subsystems 
put  on  test. 

Let  N be  the  total  number  of  a particular  noncritical  subsystem 
put  on  test.  Let  n of  these  subsystems  fail  and  their  failure  times  be 
observed.  The  remaining  m * N - n subsystems  are  removed  from  test  at 
the  truncation  points  T^,  T2,  T3»’**Tm*  Then  the  logarithmic  likelihood 
function  lnL(a,X),  based  on  the  above  sample  where  (1)  is  the  applicable 
failure  density  function,  is  given  by 


n n m 

lnL(a,X)  = nlna  + nlnX  + (a-1)  E lny.  - X E y . - X E x a (2) 

i=l  1 i=l  j=l  3 


where 


y^  •»  an  observed  failure  time 


and 


Xj  = a system  truncation  point. 


Then 


n 


N 


lnL(a,X)  « nlna  + nlnX  + (a-1)  E lnt.-X  E t° 

i=l  i=l  1 


(3) 
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» 


where 


i 


< n 


= n < i <.  N. 


This  function  yields  the  following  likelihood  equations 


and 


31nL 

3a 


n N 

- + E lnt.-X  E t“  lnt,  = 0 , 
° i=l  i-1 


(4) 


31nL  _ ii 
3X  \ 


N 

E t“  = 0 . 
i=l 


(5) 


Solving  (5)  in  terms  of  \ and  substituting  into  (4)  yields 


n 

a 


n 

+ E lnt. 
i=l 


n 


N 

E t 
i=l 


a 

i 


E t?  lnt.  = 0 . 
i=l 


(6) 


Equation  (6)  can  be  solved  for  by  a standard  iterative  procedure  but  by 
writing  (6)  in  the  form 


K(a)  = 


N 

E 

i-1 


tilnti 

"n 

i=l 


a 


n 

E lnt, 
i=l 


(7) 


8 


3 

it  has  been  shown  that  if  we  find  a K(a^)  < 0 and  a KCo^)  > 0 and  ct^ 

A 

and  a0  are  within  a sufficiently  narrow  interval  such  that  a.  < a < a„ 

a linear  interpolation  will  yield  the  required  value.  While  this  method 
will  eventually  yield  an  answer,  it  would  appear  that  the  rate  of  con- 
vergence can  be  immensely  improved  by  utilizing  basic  search  techniques. 

In  addition,  the  search  space  can  be  reduced  by  choosing  an  initial 

* 4 

that  is  relatively  close  to  a.  Dubey  has  suggested  that  the  moment  es- 

3 

timators  be  utilized  for  the  initial  value  and  Cohen  developed  a table 
of  values  relating  cx  and  the  sample  coefficient  of  variation.  Tables 
relating  the  coefficient  of  variation  to  a require  external  storage 
areas  (disks/tapes/drums)  or,  if  internal  to  a computer  program,  take 
up  precious  core  storage.  Therefore,  the  following  equation  is  utilized 
which  yields  an  initial  estimator  based  on  the  coefficient  of  variation 
which  is  within  0.1%  of  similar  table  values. 


*1  = .64364  + .18035  csch(z)  - .0317523a  + ,000684128z2  (8) 

- .00129259  csch2 (z)  + .619344(e"z)  + .534717 (e_z)4 

+ .54047 (e"Z)7 


where 


2 

g 

z = = cv  = coefficient  of  variation, 

x 

With  this  initial  estimator,  a "direct  search  with  acceleration" 
to  obtain  the  maximum  likelihood  estimate  (MLE)  of  a or  a. 


is  utilized 


3 

Cohen,  A.  C. , Maximum  Likelihood  Estimation  in  the  Weibull  Distribution 
Based  on  Complete  and  on  Censored  Samples.  Technometrics  7 (4),  579-588 
(1965). 

4 

Dubey,  S.  A.,  On  Some  Statistical  Inferences  for  Weibull  Laws,  Naval  Re- 
search Logistics  Quarterly  13,  227-251  (1966). 
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5 ct 

It  has  been  shown  that  — is  distributed  independently  of  a and  A 

/V 

and  has  the  same  distribution  as  a.  Therefore,  the  bias  in  a is  independent 
of  the  true  value  of  a and  X and  is  only  dependent  on  the  sample  size 
n.  A table  of  unbiasing  factors  B(n)  was  developed  so  that  E[B(n)a]  = a. 

The  values  in  this  table  are  approximated  by  the  following  equation: 

B (N)  = .06541717052  - .00006172777867n  + .997393979  tanh10  (n) 

-220.5624312  tanh7  (In  n/n)-1.86171021(ln  n)(tanh(ln  n/n2))  (9) 

+.01424021769 (In  n)  tanh10  (n)-39302536740000.tanh8(ln  n/n3) 

Once  the  MLE  of  a is  obtained  from  the  search  of  equation  (7),  the  ap- 
propriate unbiasing  factor  equation  (9)  is  applied  before  solving  for 
X,  using  equation  (5). 

AN  ILLUSTRATIVE  EXAMPLE 

Ten  systems  are  placed  on  test  and  the  truncation  point,  T^,  for 
each  system  is  given  (in  rounds)  below. 


T1 

- 35,142 

T6 

- 76,958 

X2 

- 43,839 

T7 

- 78,962 

T3 

- 45,540 

T8 

- 79,864 

T4 

- 61,471 

T9 

91,190 

T5 

- 66,681 

T10 

95,202 

During  the  course  of  the  test,  23  failures/replacement  (Y^)  of  sub- 
system z occurred  at  the  following  time: 

3Thoman,  D.  R. , Bain,  L.  J.,  and  Antle,  C.  E. , Inferences  on  the 
Parameters  of  the  Weibull  Distribution.  Technometrics  11  (3),  445-460 
(1969). 
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y 

1 

33,748 

y 9 

- 36,102 

yi7  - 

14,326 

y 

2 

24,898 

yio 

- 6,822 

yis  - 

27,959 

y 

3 

10,376 

yn 

- 28,824 

*19  - 

10,128 

y 

4 

27,186 

yi2 

- 37,452 

y2o  - 

31,864 

y 

5 

34,377 

yi3 

- 1,170 

*21  - 

9,366 

y 

6 

13,591 

yi4 

- 39,886 

y22  “ 

21,978 

y 

7 

878 

yis 

- 10,212 

y23  ■ 

12,924 

y 

8 

41,524 

yi6 

- 24,029 

and 

at  the 

: time 

of  each 

system  truncation,  the 

following  subsystem  times 

(xp 

were 

observed: 

X1  " 

1,395 

x 6 

40,855 

X 2 " 

8,564 

X 7 ‘ 

5,863 

X 3 ' 

18,354 

X 

8 

38,806 

X 4 ' 

12,624 

X 9 " 

14,662 

x - 

5 

25,157 

x ■ 

10 

8,940 

This  sample  for  subsystem  z is  from  a population  in  which  a = 1.364  and 

X = .00000102336.  This  data  can  be  summarized  as: 

23  23  . 

n = 23 , N=33,  Z t. =499618,  Z t =14396873864. 

1 1 

It  follows  that  x = 21722.52,  s^=  154083086.858,  and  the  coefficient  of  variation 
cv  * .3265.  Then  from  (8)  the  initial  estimate  of  a=1.81.  To  obtain  the 
maximum  likelihood  estimate,  a is  varied  in  equation  (7)  until  |K(a)|  is  mini- 
mized.  This  results  in  a=1.6421  and  applying  the  unbiasing  factor  for  a 

/v 

sample  size  of  33,  a = 1.5731. 
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TEST  OF  HYPOTHESIS 


Consider  the  following  failure  density  function: 

ft— 1 —X 

f(t;a,X)  * aXt  e dt  a>0;X>0;t>0  . 

When  the  resulting  sample  likelihood  function  of  a life  test  can  be 

then,  this  function  yields  the  following  partial  derivatives  of  the 
rithmic  likelihood  equations: 


91nL 

9a 


n N 

— + £ lnt.  - X I t “ lnt.  = 0 , 

° i-1  1 i-1  1 1 


91nL  _ n 
9X  X 


- 0 , 


92lnL 


9a 


~ X E t.°  (lnt.)2  = 
a i=l 


0 , 


92lnL 


9X 


0 , 


91nL 

9a9X 


N 

I t.°  lnt.  = 0 • 
i=l 


(10) 

written 

(ID 

loga- 

(12) 

(13) 

(14) 

(15) 

(16) 
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The  maximum  likelihood  estimate  (MLE)  a and  X can  be  obtained  by 
solving  (13)  in  terms  of  X and  substituting  into  (12) . Having  obtained 
the  MLE  for  the  parameters  of  the  Weibull  law,  it  is  then  possible  to 
determine  confidence  limits  for  meaningful  parametric  functions  such  as 
durability  and  reliability. 

Durability  can  be  defined  as  the  probability  that  a randomly  selected 
item  from  an  infinite  lot  will  continue  to  perform  satisfactorily  without 
a durability  failure  beyond  to  and  is  given  by 

oo 

Prob[T  _>  to]  * / f1p(t)dt  = exp(-Xt0a)  • (17) 

to 

If  the  parameters  a and  X are  known,  the  above  problem  is  completely  solved 
and  the  exact  answer  is  given  by  (17);  however,  usually  these  parameters 
are  unknown.  It  is  well  known  that  the  large  sample  MLE  is  approximately 
normally  distributed  about  the  true  parameter  value  as  a mean  for  large 
samples.  This  is  a powerful  tool  and  will  be  used  to  establish  confidence 
limits  on  durability  when  the  true  parameters  a and  X are  unknown.  Assum- 
ing the  above,  it  follows  that^ 


Ci  2 
C22 


a-a 

A 

x-x 


X2(2) 


(18) 


where 


Cn 


(19) 


Mood,  A.  M.,  Introduction  to  the  Theory  of  Statistics.  McGraw-Hill  Book  Co., 
Inc.,  New  York,  New  York,  1950. 
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(20) 


Czz 


■ -Kf]  • 


and 


From  (14)  and  (19) 


Cl2  “ C21 


_E  f 3 2 Ink] 

E [_  3a3xJ 


Cl  1 - E 


[?  + \-i  ti,’<lnti)2]  • 


(21) 


(22) 


Examine 


f t°(lnt) 

J 0 


2 . a-l-Xt  , 
aXt  e dt 


and  let 


b) 

db) 

lnu 

lnt 

2 


(lnt) 


t° , 

ata-1dt 


alnt  , 

- lnu  , 
a * 

— —(lno))^ 
a 


then 


k 2 — Xu>  , 

e da) 


I * -4r  f (DO-nu))' 

a2  -A) 

From  page  578  of  referenced  literature7,  equation  4.358(2) 


7Gradshteyn,  I.S.,  and  Ryzhik,  I.M.,  Table  of  Integrals,  Series,  and  Products. 
Academic  Press,  New  York,  New  York.  1965. 


14 


r xv  Xe^x (lnx) 2dx 


r(v)/uV{[iKv)-ln(u)]2+  5 (2,v-l) } 


y > 0,v  > 0 


now  let 

x * w, 

V - 2, 

and 

y * X, 

then 

X z'  oo 

1 " ~J  0)(lna))2e~Xwdw  - Xr(2)/ (a2X2) { [\|/(2)-lnX]2+^(2 ,1) } 

a 

= {[iK2)+ln  -y-]2  + 5(2,1)} 

Xot 


Q 

From  page  1073  of  referenced  literature  , equation  9.521(1)  the  Riemann's 
Zeta  Function 


5(z,q)  = E — j z > 1 

N=0  (q+N) 


•*.  5(2,1) 


g 

From  page  260  of  referenced  literature  , equation  6 •4*10  the  Polygamma 
Function 


i|)N(z)  = (-1)N+1N!  E (z+k)  N_1  (z^0,-l,-2, . . . ) 

k=0 


g 

Abramowitz,  M. , and  Stegun,  I.  A.,  Handbook  of  Mathematical  Functions  with 
Formulas,  Graphs,  and  Mathematical  Tables,  National  Bureau  of  Standards.  US 
Government  Printing  Office,  Washington,  DC.  1970 
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/(l)  - (-1)  1!  I (l+k)-2=  j + + “V  + ~T  + 

k«0  1 2Z  T S 


• • • ) 


therefore. 


C(2 ,1)  « r(l) 


and 


1=  {[<K2)+ln  ^]2  + /(l)}  . 

Xa2  A 


Equation  (22)  can  now  be  written  as 


C11  = "T  + "T  {W<2>+ln  X ]2  + ^(D) 
a a 


(23) 


where 


and 


i|»  (2)  = .4227843351  , 


\p  (1)  = 1.6449340668 


From  (16)  and  (21) 


Examine 


:12  ' {f=1  '>‘i] 


r a , ,„cx-l  -Xt“ 

/ t lnt  aXt  e dt 

Jo 


(24) 


and  let 


0)  « 

rt 

P 

do)  = 

at0'1, 

lnco  = 

alnt 

lnt  = 

lnw/a 

16 


then 


- r 

a J o 


wlnu)  e 


-Aw 


From  page  576  of  referenced  literature  , 


du  . 

equation  4.352(1) 


f xv  1 e PX  lnxdx  = — r(v)  [if»(v)-ln(y)  ] u > 0,v  > 0 
JO  „v 


let 


then 


v = 2 , 

P = A , 

x = (1)  , 


I - £ /*“  wlnw  e Au)  du  = (£)(-4)r(2)[iK2)-lnXJ 

a Jq  a 

■ dr b ■ 


Equation  (24)  can  now  be  written 

C12  ' h i]  , 

and  it  follows  that 


(25) 


(26) 


As  just  shown, 


the  C 


i.j 


matrix  is  the  inverse  asymptotic  variance- 


covariance  matrix  of  (A, a)  and  is  obtained  by  taking  the  negatives  of 
the  expected  values  of  the  second  order  derivatives  of  logarithms  of  the 
likelihood  functions.  Using  (18),  we  can  now  obtain  the  appropriate  con- 
fidence limits  for  the  true  probabilities  . In  the  case  of  reliability 


Loc. 


Cit. 
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or  durability,  we  want  to  insure  that  the  true  reliability  or  durability 

is  above  some  minimum  value  R or  D at  a desired  confidence  level.  From 

m m 

(18),  then, 

Cn(a-a)2  + 2C12(X-X)  (a-a)  + (X-X)=y2(2)  . (2?) 


The  above  equation  describes  the  boundary  of  a confidence  region  in  the 
parametric  space  (a,X).  It  is  approximately  chi-squared  distributed  with 
two  degrees  of  freedom  and  will  serve  to  determine  an  ellipsoidal  confi- 
dence region  in  the  (a,X)  space.  In  order  to  obtain  the  upper  confidence 
limit,  all  that  is  required  is  to  search  the  parameter  space,  (27), 
varying  a and  X until  the  durability  function  (17)  is  maximized. 


A PROCEDURE  FOR  COMPUTING  EXPECTED  NUMBER  OF  RENEWALS 

Let  X^,  X2,  X^,...  be  a renewal  process,  that  is,  a sequence  of 
independent,  nonnegative  and  identically  distributed  random  variables 
which  are  not  all  zero  with  probability  one,  with  the  probability  den- 
sity function  f(X)  and  the  distribution  function  F(X).  = X^+  X2+. . . 

X^  is  interpreted  in  renewal  theory  as  the  time  up  to  the  kth  renewal 
and  the  probability  that  < x is  given  by  the  k-fold  convolution  of 


F(x) 


(x-t)  d F(t) , 


where 

Fq (x)  = 1. 

The  primary  purpose  of  this  procedure  is  to  calculate  N^  which  is 
defined  to  be  the  maximum  suffix  K for  which  t,  subject  to  the 


18 


convention  N » 0,  if  > t.  In  this  application,  the  X.^  represents 
successive  lifetimes  of  the  object  being  renewed,  and  N^  is  the  number  of 
renewals  made  by  time  t,  subject  to  the  original  object  having  been  in- 
stalled at  time  0,  i.e.,  N£  is  the  number  of  renewals  in  (0,t). 

The  usual  procedure  for  determining  various  functions  of  Renewal 
Theory  is  to  find  the  Laplace  Transforms  of  these  functions  and  then  re- 
vert to  the  time  domain.  However,  in  the  case  of  the  Weibull  distribution, 

9 

this  approach  is  not  convenient  and  Smith  and  Leadbetter  have  developed 
an  expansion  of  the  Renewal  Function  into  a power  series  of  tawhere  a is 
the  Weibull  shape  parameter.  White^  demonstrated  a procedure  for  eval- 
uating the  higher  moments  and  cumulants  of  the  number  of  renewals  Nfc  and 

1/2 

computed  the  mean  M(t)  and  the  standard  deviation,  (Var  N^.)  , of  N^ . 

2 

Lomnicki  developed  a method  to  evaluate  the  distribution  of  the  number  of 
renewals  which  is  defined  by  the  family  of  functions  Wk(t)  where 

Wk(t)  = Fk(t)  ~ Fk+l(t)  (k=0 , 1 , . . . ) 

and  is  the  probability  of  exactly  k renewals  in  (0,t). 

In  order  to  be  able  to  evaluate  more  of  the  various  forms  of  the 
renewal  functions  in  the  case  of  the  Weibull  renewal  process,  the  basic 
method  developed  by  Lomnicki  was  used. 

As  shown  by  Lomnicki,  if  we  assume  that  Wk(t)  is  represented  by  the 
unique  series 

<30 

Wfc(t)  = E ak(s)Ps(ta) 
s=k 

9 

Smith,  W.  L. , and  Leadbetter,  M.  R. , On  the  Renewal  Function  for  the 
Weibull  Distribution.  Technometrics  5(3),  393-396  (1963). 

F^White,  J.  S.,  Weibull  Renewal  Analysis.  Third  Annual  Aerospace  Reliability 
and  Maintainability  Conference,  pp  639-65/ . Society  ot  Automotive 
Engineers,  Inc.  New  York,  New  York.  1964. 

2 

Loc.  Cit. 
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where 


then 


pk(o 


(k-0,1,2, . . .)  , 


Vt) 


Z 

8“k 


Ps(ta) 


E (-l)P+k0 

p=k  p 


bk(p) 

y(p) 


9 


and  we  have 


Vs) 


i (-i)p+lt 

P“k 


bk(P) 


(S) 

V Y (p) 


(k=0,l, . . . ;s=k,k+l, . . . ) 


where 


y(r)  * r (ar+1) /T (r+1)  (r-0,1,...) 


and 

bg(s)  = Y(s) . (6=U , 1 , . . . ) 


s-1 

b(s)  = Z b.  (r)Y(s-r)  (k=0,l, . . . ;s=k,k+l,k+2, . . .)  . 

k+1  r=k  k 


The  expected  number  of  renewals  M(t)  = E[N(t)]  which  is  the  renewal 
function  is  given  by 


Since 


is  equal  to 


OO 

M(t)  = Z n P[N(t)=n] 
n=0 


00 

M(t)  = Z n W (t) 
n n 


Wx(t)  + 2W2(t)  + 3W3(t)  + ... 


F1(t)  -F2(t)  +2F2(t)  -2F3(t)  + 3F3(t)  -3F4(t)  + ...  , 
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then 


CO  00 

l n W (t)  - Z F (t) 
n=0  n=0 

and 

00 

M(t)  = Z Fn(t)  , 


where  Fn(t)  is  the  convolution  of  the  cumulative  distribution  function 
F(t)  of  the  first  renewal  time  and  the  n-1  subsequent  renewal-time  dis- 


tributions and  is  the  distribution  of  S , the  time  of  the  nth  renewal. 

n’ 

Now, 


Fk(t)  - Hk(t>  + Fk+1(t) 

- ...  + Wk+J(t)  + Fk+3(t) 

■ 

• » 

Fk(t)  - Z Wr(t)  (k=l,2,3. . .) 

r=k 


or 

k-1 

Ffc(t)  = 1 - Z W (t), 
r=0 


so  we  can  now  express  Fk(t)  as 


Fk(t)  = Z Z ar(s)  P8(ta) 
r=k  s=r 


Substituting  the  Poisson  cumulative  function 


Ds(t)  - 

r=s 
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means  that 


00 


where 


and 


so  that 


or 


where 


Fk(t)  = Z 6k(s)  Ds(t°) 
s“k 


pk(k)  = \(k)  ’ 


^(s)  = Z ar(s)  - Z a (s-1)  (s  > k)  , 
r=k  r=k 


oo  oo  s 

H(t)  = Z F (t)  *=  Z D (t°)  Z B (s)  , 
k=l  K s=l  8 k»l  k 


M(t)  = Z C(s)  D (ta)  , 
k=l  8 


C(s)  = Z g (s)  . (s=l,2,...)  • 

k=l  k 


We  can  now  evaluate  M^(t),  the  Renewal  Function,  W^(t);  the  probability 
of  exactly  k renewals;  and  F^(t),  the  distribution  of  time  of  the  kth 
renewal.  These  functions  are  very  useful  in  evaluating  models  of  relia- 
bility, inventory,  and  queueing  process. 

V‘>  * E<T>  ■ y£>  <Fn(t,-Fn+l(t»  ' 

n=u 

The  moments  of  N(t)  may  be  derived  from  M^Ct)  by  the  relationship 

k k 

E (N (t ) ) - Z t.  M (t)n! 
n-1  k’n  n 
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where  t^  n are  the  Stirling  numbers  of  the  second  kind  and  represent  the 
number  of  ways  of  partitioning  a set  of  k elements  into  n nonempty  subsets. 
In  particular 

E{N(t)  } = M^t)  , 

and  9 

E{N(t)  j - M^(t)  + 2M2(t)  . 

Or,  if  desired,  the  moments  of  N(t)  may  be  derived  from  the  definition  of 
expected  values. 

GO 

E{N(t)K}  ="  I nK  Prob(N(t)=n)  , 
n*0 


E{N(t)K}  ■=  Z 

"k<f„  • 

n*0 

00 

E{N(t)K}  = Z 

nSf  (t)  , 

n=0 

n 

00 

E{N(t)K}  = Z 

nK-(n-l)KFn(t)  . 

Next  page  is  blank. 


23 


LITERATURE  CITED 


1.  Mazza,  Thomas  M.  and  Banash,  Robert  C.,  PAA-TR-1-73. 

Decision  Risk  Analysis  for  XM204,  105MM  Howitzer,  Toved  Reliability/ 
Durability  Requirements,  Systems  Analysis  Division,  Plans  and  Anal- 
ysis Directorate,  US  Army  Weapons  Command,  April  1973. 

2.  Lomnicki,  Z.  A.,  A Note  on  the  Weibull  Process.  Biometrika 
53,  375-381  (1966). 

3.  Cohen,  A.  C.,  Maximum  Likelihood  Estimation  in  the  Weibull 
Distribution  Based  on  Complete  and  on  Censored  Samples.  Technometrics 
7 (4),  579-588  (1965). 

4.  Dubey,  S.  A.,  On  Some  Statistical  Inferences  for  Weibull 
Laws,  Naval  Research  Logistics  Quarterly  13,  227-251  (1966). 

5.  Thoman,  D.  R.,  Bain,  L.  J. , and  Antle,  C.  E. , Inferences  on 
the  Parameters  of  the  Weibull  Distribution.  Technometrics  11  (3) , 
445-460  (1969). 

6.  Mood,  A.  M. , Introduction  to  the  Theory  of  Statistics. 
McGraw-Hill  Book  Co. , Inc.  New  York,  New  York,  1950. 

7.  Gradshteyn,  I.  S.,  and  Ryzhik,  I.  M. , Table  of  Integrals, 
Series,  and  Products.  Academic  Press,  New  York,  New  York.  1965. 

8.  Abramowitz,  M. , and  Stegun,  I.  A.,  Handbook  of  Mathematical 
Functions  with  Formulas,  Graphs,  and  Mathematical  Tables,  National 
Bureau  of  Standards.  US  Government  Printing  Office,  Washington,  DC. 
1970. 


9.  Smith,  W.  L.,  and  Leadbetter,  M.  R. , On  the  Renewal  Function 
for  the  Weibull  Distribution.  Technometrics  5 (3) , 393-396  (1963) . 

10.  White,  J.  S.,  Weibull  Renewal  Analysis.  Third  Annual  Aero- 
space Reliability  and  Maintainability  Conference,  pp  639-657.  Society 
of  Automotive  Engineers,  Inc.  New  York,  New  York.  1964. 


Next  page  is  blank. 


25 


DISTRIBUTION  LIST 


No.  of  Copies 
12 


1 


1 


1 

1 

1 

1 


1 


1 


Organization 


Commander 

Defense  Documentation  Center 
ATTN:  TIPCR 

Cameron  Station 
Alexandria,  VA  22314 

Commander 

US  Army  Materiel  Command 
ATTN:  AMCDL 

5001  Eisenhower  Avenue 
Alexandria,  VA  22304 

Commander 

US  Army  Materiel  Command 
ATTN:  AMCMA 

5001  Eisenhower  Avenue 
Alexandria,  VA  22304 

Commander 

US  Army  Materiel  Command 
ATTN:  AMCRD 

AMCRD-R 

AMCRD-U 

AMCRD-W 

5001  Eisenhower  Avenue 
Alexandria,  VA  22304 

Commander 

US  Army  Materiel  Command 
ATTN:  AMCQA 

5001  Eisenhower  Avenue 
Alexandria,  VA  22304 

Commander 

US  Army  Materiel  Command 
ATTN:  AMCPA-SA 

5001  Eisenhower  Avenue 
Alexandria,  VA  22304 


27 


No.  of  Copies 


Organization 


1 

2 

1 

1 

1 

1 


1 

1 

1 


1 


1 


1 


1 


1 


1 


Commander 

US  Army  Armament  Command 
ATTN:  AMSAR-MAE 

AMSAR-SAS 
AMSAR-RD 
AMSAR-QAR 
AMCPM-CAWS 
AMSAR-CP 

Rock  Island,  IL  61201 

Director,  Army  Materiel  Systems  Analysis  Agency 
Aberdeen  Research  and  Development  Center 
ATTN:  AMXSY-E 

AMXSY-M 
AMXSY-R 

Aberdeen  Proving  Ground,  MD  21005 
Commander 

US  Army  Aviation  Systems  Command 
ATTN:  AMSAV-D 

PO  Box  209 

St.  Louis,  MO  63166 
Commander 

US  Army  Ballistics  Research  Laboratory 
Aberdeen  Proving  Ground,  MD  21005 

Commander 

Harry  Diamond  Laboratories 
ATTN:  Systems  Analysis  Office 

Washington,  DC  20438 

Commander 

US  Army  Electronics  Command 

ATTN:  AMSEL-SA 

Fort  Monmouth,  NJ  07703 

Commander 

US  Army  Electronics  Command 

ATTN:  AMSEL-PL 

Fort  Monmouth,  NJ  07703 

Commander 

US  Army  Mobility  Command 

ATTN:  AMSME-XS 

4300  Goodfellow  Blvd. 

St.  Louis,  M0  63120 


28 


# 


No.  of  Copies 


Organization 


1 


1 

1 


1 

1 

1 


1 


1 

1 

1 

1 


1 


1 

1 

1 


Commander 

US  Army  Missile  Command 

ATTN:  AMSMI-DA 

Redstone  Arsenal,  AL  35869 

Commander 

US  Army  Tank-Automotive  Command 
ATTN:  AMSTA-RPL  (Tech  Lib) 

AMSTA-SR 

Warren,  MI  48090 
Commander 

US  Army  Test  and  Evaluation  Command 
ATTN:  AMSTE-SY 

AMSTE-RA 

Aberdeen  Proving  Ground,  MD  21005 
Commander 

US  Army  Troop  Support  Command 
ATTN:  Systems  Analysis  Office 

4300  Goodfellow  Blvd 
St.  Louis,  MO  63120 

Commander 

US  Army  Safeguard  Evaluation  Agency 
ATTN:  SSEA-T 

White  Sands  Missile  Range,  NM  88002 
Commander 

US  Army  Training  and  Doctrine  Command 
ATTN:  ATCD-CI 

ATCD-CF 
ATCD-CP 
ATCD-SI 

Fort  Belvoir , VA  22060 
Commander 

USA  Combined  Arms  Combat  Developments  Activity 
ATTN:  ATCACC 

Fort  Leavenworth,  KS  66027 
Commander 

Rock  Island  Arsenal 
Rodman  Laboratories 
ATTN:  SARRI-LA 

SARRI-LR 
SARRI-LE 

Rock  Island,  IL  61201 


29 


No.  of  Copies 


Organization 


1 


1 


5 


2 


2 


1 


* 


Commander 

Watervliet  Arsenal 
ATTN:  SWEWV-RDD 

Watervliet,  NY  12189 

Manager,  Defense  Logistics  Information  Exchange 
USA  Logistics  Management  Center 
ATTN:  Mr.  Leon  J.  Scarbrough 

Fort  Lee,  VA  23801 

Director 

US  Army  Management  Engineering  Training  Agency 
ATTN:  AMXOM-AMS 

Library 

Rock  Island,  IL  61201 
Commander 

US  Army  Logistics  Management  Center 
ATTN:  AMXMC-LS 

Fort  Lee,  VA  23801 

Commander 

US  AMC  Intern  Training  Center  - USALMC 
ATTN:  AMXRR-IT 

Texarkana,  TX  75501 

Commander 

Rock  Island  Arsenal 
ATTN:  SARRI-LPL 

Rock  Island,  IL  61201 


i 


30 


- 


